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ABSTRACT 

Air-ffuorescence detectors such as the High Resolution Fly's Eye (HiRes) 
detector are very sensitive to upward-going, Earth-skimming ultrahigh energy 
electron-neutrino-induced showers. This is due to the relatively large interaction 
cross sections of these high-energy neutrinos and to the Landau-Pomeranchuk- 
Migdal (LPM) effect. The LPM effect causes a significant decrease in the 
cross sections for bremsstrahlung and pair production, allowing charged- current 
electron-neutrino-induced showers occurring deep in the Earth's crust to be de- 
tectable as they exit the Earth into the atmosphere. A search for upward- going 
neutrino-induced showers in the HiRes-II monocular dataset has yielded a null 
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result. From an LPM calculation of the energy spectrum of charged particles as 
a function of primary energy and depth for electron-induced showers in rock, we 
calculate the shape of the resulting profile of these showers in air. We describe 
a full detector Monte Carlo simulation to determine the detector response to 
upward-going electron-neutrino-induced cascades and present an upper limit on 
the flux of electron-neutrinos. 

Subject headings: cosmic rays — neutrinos — acceleration of particles — large- 
scale structure of the universe 



1. Introduction 



We report on a search for upward-going electron-neutrino showers in the High- Resolution 
Fly's Eye II data set, and on the upper limit on the flux of u e set by the HiRes-II detector . 
The HiRes project has been discussed previously (lAbu-Zayyad et al.lll999uBoyer et al.ll2002l ); 
the detector is an air-fluorescence detector located on two sites 12.6 km apart in Utah at 
the U.S. Army Dugway Proving Ground. The HiRes-II detector, located on Camel's Back 
Ridge, is composed of 42 spherical mirrors of 3.7 m 2 effective area covering nearly 360° in 
azimuth. Half of these, known as ring-one mirrors, cover between 3°-17° in elevation; the 
other half (ring-two) cover between 17°-31° in elevation. 

Cosmogenic neutrinos, with energies mostly in excess of 10 18 eV, are produced via tt 
and \i decays following photopion production fro m high-energy cosmic ray protons incident 
on the cosmic microwave background radiation (jSteckerlll968t iMargolis et al.lll978l ). There 
is evidence to suggest that gamma-ray bursts an d active galactic nuclei je t s are possible 
sourc es of high-energy cosmic rays and neutrinos (IWaxman fc Bahcalll 119971 ; lHalzen fc Zas 
1997). Several theoretical limits on the flux of cosmogenic neutrinos have been proposed 
Semikoz fc Sirflbooi ISeckel fc Standi l2005h . 



Although large uncertainties exist, neutrino cross secti ons have be en calculated to vary 
from ~ 1(T 32 cm 2 at 10 18 eV to ~ 1(T 31 cm 2 at 10 21 eV (lRenoll2005h . The opacity of the 
earth to neutrinos at these high energies therefore prohibits the detection of any upward- 
going event with an elevation angle larger than a few degrees. 

In the charged- current interaction of a v e in the earth's crust, a high-energy electron will 
be created. The electromagnetic cascade generated by the electron will develop much more 
slowly due to the onset of the Landau- Pomeranchuk-M i gdal ( LPM) effect. The LPM effect, 
first described classically by lLandau fc Pomeranchukl (119531 ) and later given a quantum- 
mechanical treatment by Migdal ( 19561 ). predicts that the cross sections for bremsstrahlung 
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and pair-production should decrease for a high-energy charged particle propagating in a 
dense medium, effectively slowing and elongating the develop ment of the result i ng sh ower 
of particles (a detailed , more modern approach can be found in iTakahashi et al.l (120031 ) and 
Baier & Katkov! (120041 ) ) . The energy at which this effect becomes appreciable is inversely 
proportional to the square of the Lorentz factor 7, and therefore the LPM effect should be 
much more pronounced for the showers generated from a v e charged- current interaction than 
for showers precipitated by or v T in the energy range in which HiRes is sensitive. 

It is most probable that a neutrino-induced electromagnetic cascade would be long and 
nearly-horizontal and observed primarily in the HiRes-II ring-one mirrors. Due to the LPM 
effect, one expects electron-neutrino- induced showers that begin several tens to hundreds of 
meters deep in the crust to emerge with enough charged particles to be detected by HiRes-II, 
thereby increasing the effective aperture of the detector at high energies. 



2. Search for upward-going neutrino events 

The entire HiRes-II data set, which extends from late 1999 to Spring 2006, was con- 
sidered when searching for evidence of neutrino-induced upward-going showers. Using the 
standard routines that were developed for analyzing downward-going cosmic-ray events, we 
reconstructed the trajec t ories of each upward-going event based on the measured timing and 



geometry (see ISokolskyl (Il989l ) for a description of time- and plane-fitting for extensive air 
showers) . 

The data were then filtered in time and position to exclude all calibration laser events, 
which resulted in a loss in the detector aperture of less than 1%. Additionally and consis- 
tent with standard procedure for the analysis of cosmic-ray data, events were rejected that 
passed within 100 meters of the detector, had track lengths smaller than 10°, and that had 
geometrical uncertainties from timing greater than 36°. 



3. The Landau-Pomeranchuk-Migdal Effect 



At electron energies below the LPM threshold energy (61.5 L cm TeV (IStanev et al.lll982l ). 
where L cm is the interaction length in cm), the longitudinal profile of an electromagnetic 
shower can be well approximated by the relation 
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This functional form was first described by iGreisenl (119561 ). with /3q as the log of the ratio 
of the energy of the incident electron to its critical energy E c , t as the depth in radiation 
lengths and s = u^Po] • 



This relation begins to break down at high energies, greatly underestimating the distance 
over which the electromagnetic cascade evolves due to the decrease in the cross sections for 
bremsstrahlung and e + e~ pair production. Studies of the electron shower profiles i n rock, 
water and lead above the LPM threshold ene rgy have been conducted previously (IMisaki 



19901 ; IStanev et al.l 11982 ; lAlvarez-Munia Il999l ) . As expected, the results of these analyses 



show that the shower profiles of electron-induced cascades are elongated significantly with 
respect to the Greisen approximation at energies above the LPM threshold, and evolve 
differently based on the densities of the media in which the showers propagate. 



4. Calculation of sensitivity to electron-neutrino showers 

To simulate z/ e -induced electromagnetic cascades, we used a four-step process. First, 
we calculated the average profiles of electron- induced showers using the LPM effect. We 
then used a Monte Carlo method to simulate the arrival directions and interaction points of 
u e around the HiRes detector. The shower profiles in air were then passed into the HiRes 
detector Monte Carlo to calculate the amount of light seen by the detector. The HiRes 
analysis programs were then run on the resulting Monte Carlo events to arrive at a v e 
aperture. 



4.1. Calculating electron- neutrino-induced electromagnetic cascade profiles 

In order to treat charged-current v e N interactions in the earth's crust, it is necessary to 
understand the physics of the transition of an electromagnetic cascade from a dense medium 
to a less dense medium (namely, from rock to air). It is therefore important not only to know 
the number of charged particles after traversal of a given amount of material in rock, but also 
the energy spectrum of these particles as they leave the ground and enter the atmosphere. 



We followed the formalism of IStanev et al.l (119821 ) for calculating the energy-dependence 



of the probabilities for undergoing pair production and bremsstrahlung at LPM energies. 
Taking into account any other losses (e.g. Compton scattering and ionization energy loss), 
we calculated two functions: N™ ck (E , E,d) and N™ r (E ,E,d), which describe the average 
number of charged particles with energy E resulting from the cascade of an electron, positron 
or photon with initial energy Eq after traversing an amount of material X in rock or air. The 
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functions N r e ock and Nf r were determined for Eq at every decade between 10 12 and 10 21 eV 
using our LPM calculation and from E c to 10 12 eV using Equation [TJ LPM calculations of 
shower profiles from particles with initial energies below 10 12 eV were found to be nearly 
identical to profiles calculated using Equation [TJ 



4.2. Simulating neutrino events 

We approximated the earth as a sphere with a radius equal to that at the Dugway 
Proving Ground in Utah. The density below 58.4 km beneath the surface (mantle) and the 
density from from 58.4 km to the surface (crust) were taken to be 4.60 and 2.80 g cm -3 
respectively. The atmosphere was also simulated up to a height of 50 km above sea level. 

Electron-neutrino energies were considered from log E v of 18 to 21. The energy-dependence 
and inelasticity of the charged- and neutra l-current vN in t eraction cross sections wer e calcu- 



lated based on the pQCD CTEQ5 model ([Lai et al.ll2000l ; iGazizov fc Kowalskill2005l ). From 



the ratio of the cross sections for charged-current (CC) and neutral-current (NG) interac- 
tions, 70% of the events were thrown as CC events, while the remaining 30% were considered 
NC events. 

Neutrino arrival directions were chosen at random such that they only penetrated the 
atmosphere no more than 15° below the horizon. Events with elevation angles greater than 
15° do not contribute appreciably to the HiRes-II total u e aperture due to the very small 
probability of their transmission through the crust and mantle and subsequent interaction 
near the detector (a 10 18 eV neutrino at 15° has a probability of ~ 10~ 12 of transmission 
and interaction near the detector; this value drops to ~ 10~ 60 at 10 21 eV). The variables 
describing the geometry of the neutrino trajectory were determined, such as the distance of 
closest approach to the detector, the vector normal to the shower-detector plane, and the 
angle of the shower in the shower-detector plane. 

For these earth-skimming events, the traversal of a critical amount of material X c (mea- 
sured in g cm -2 ) was found such that when the shower emerges from the rock into air, it 
contains at least 10 7 charged particles; showers with a maximum number of charged particles 
less than 10 7 will not trigger the HiRes detector. This critical pathlength is used to separate 
the probabilities for neutrino transmission and interaction. The transmission probability e t , 
was calculated as the probability for a neutrino to penetrate up to X c . The interaction prob- 
ability Si, was calculated from the pathlength of the neutrino from X c until escape from the 
atmosphere. Since the amount of material traversed in the interaction region is always much 
less than the mean neutrino interaction length, the actual point of interaction for each neu- 
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trino was then chosen at random for distances X > X c . For neutrinos with small elevation 
angles that do not pass through the Earth, we considered events that entered the atmosphere 
above the horizon as well as those that interacted below the horizon and yielded at least 
10 7 particles at the horizon. For events interacting below the horizon, X c was taken to be 
amount of air penetrated at the horizon. In the case of events that entered the atmosphere 
above the horizon, we set e t to unity and calculated from the total distance traversed in 
the atmosphere. 

For all z/ e N interactions the energy transferred to the secondary electron or hadron was 
chosen from the inelasticity distribution (da/dy) for the pQCD CTEQ5 model. For Earth- 
skimming CC events, the resulting observable profile in air was found from a superposition 
of showers obtained from the energy spectrum of electrons, positrons and photons emerging 
from the rock. The profiles of CC events that did not pass through the earth were interpolated 
from the N™ r functions described in the previous section . The profile s for all NC events 



were calculated using the standard Gaisser-Hillas model (jGaisserl Il990l ). Each profile was 
then weighted by a factor w = e t £i, to describe the total probability of transmission and 
interaction near HiRes-II. Figure Q] shows the average profiles of five electron-induced air 
showers emerging from the ground at different depths along an average 10 20 eV electron- 
induced shower in rock. 



4.3. Simulating detection by HiRes 

Having generated shower profiles using the LPM effect, the shower profiles were then 
passed through a HiRes Monte Carlo program which models the response of the detector 
to cosmic-ray-induced showers. This program determines the amount of fluorescence and 
Cerenkov photons produced for a given number of charged particles, and scatters and at- 
tenuates the light appropriately when given the known variables describing the geometry 
of the shower with respect to the detector. The program then models the HiRes-II trigge r 



conditions to decide if the simulated shower is read out by the detector (lAbbasi et al. 112004 ). 



4.4. Analysis and filtering of simulated events 

Simulated showers which triggered the detector were analyzed using the same routines 
used in the analysis of the real data used in our search for neutrinos in the upward-going 
HiRes-II data. The variables describing the geometry of the shower were fit and compared 
to the known variables. An event was considered accepted when it passed the same cuts 
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described in Section 2. 



5. Calculating an electron-neutrino flux upper limit 

For the purposes of arriving at a predicted HiRes-II u e aperture, the simulated events 
were collected in 30 0.1-decade energy bins from 10 18 to 10 21 eV. The aperture for a given 
energy bin was found 
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where R is the radius of the earth extended 50 km to the edge of the atmosphere. The 
geometrical component of the aperture is derived from the area and solid angle of a 30° cap 
on a sphere of radius R. This is then adjusted by the weighting factor w (discussed in Section 
4.2) for each of the Na events that trigger the detector out of a total Nt events thrown in 
the given energy bin. The HiRes-II v e aperture is shown in Figure [2J 



Consistent with our study of v T ( ]Martensll2007l ). we calculate a flux limit in three energy 
bins: AE = 10 18 - 10 19 , 10 19 - 10 20 , and 10 20 - 10 21 eV, over the total HiRes lifetime of 
3638 hours. We observe no neutrino events over the entire energy range. We calculate the 
flux limit (E 2 ^) at the 90% confidence level to be 4.06 x 10 3 , 3.55 x 10 3 and 4.86 x 10 3 eV 
cm -2 sr _1 s" 1 at 10 18,5 , 10 19 ' 5 and 10 20 ' 5 eV, respectively. Combined with our v T results and 
provided equal mixing of all neutrino flavors, this reduces the limit to 3.81 x 10 2 , 9.73 x 10 3 
and 4.71 x 10 3 eV cm -2 sr" 



1 s" 1 . 



6. Discussion 

As is the case with all high-energy neutrino calculations, the largest uncertainty lies in 
the extrapolation of vN cross sections. Different cross section models can cause the limits to 
vary somewhat. The incorporation of cross sections from previous and more recent versions 
of the CTEQ model can change the limits by as much as 10 to 40% at the lowest and highest 
energies, respectively. 

Recent work imposing the Froissart bound on structure functions for extrapolating vN 
cross sections show a decr ease in cross sections at 10 21 eV by about a factor of 8 over the 
CTEQ5 parameterization (jBlock! 120071 ). These cross sections increase our v e limit by 40% 



at the lowest energy bin and increase the value of our highest energy bin by a factor of 3. 
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In addition to uncertainties in vN cross sections, our limits are also sensitive to the 
energy transferred to the second ary electron or had ron. From parameterizations of the mean 
inelasticity in vN interactions (jOuigg et al.lll986l ). if we allow the transfer of exactly 80% 
of the neutrino energy to the electron (and 20% to the hadron), our v e limits will increase 
between about 15% at 10 18 - 5 eV, remain unchanged in the middle energy bin and decrease 
by about 5% at 10 21 - 5 eV. 



7. Conclusion 



We have found no evidence of upward-going neutrino-induced cosmic-ray showers in 
the HiRes-II data. We have presented a technique for modeling the full HiRes-II detector 
response to ultrahigh energy neutrino-induced LPM cascades in rock and air. With no 
neutrino events seen in the HiRes-II data, and provided equal mixing of all neutrino flavors, 
we have found an upper limit on the flux of ultrahigh energy neutrinos at a 90% confidence 
limit. 

Figure [3] shows the upper limit on the neutrino flux from the analysis of the HiRes v e 
and v r flux limits as compared to three theoretical curves and to calculated flux limits from 
other experiments. The v e flux limits reported here have improved upon those for the Fly's 
Eye by about two and a half orders of magnitude. Combined with the results of t he v T 
analysis, this limit lies just above the theor etical neutrino flux of ISemikoz fc Sigll (120041 ). and 
about an order of magnitude above that of ISeckel &: Stanevl (120051 ). Our combined neutrino 
flux limit is a bout two and a half orders of magnitude above the cosmogenic neutrino flux 
predictions of iBrusoval (120071 ). which has been derived from a proton injection model with 
cosmologically evolving sources and injection spectra that fit the HiRes cosmic-ray spectrum. 
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Fig. 1. — An average 10 20 eV electron shower profile in rock (solid line) with average shower 
profiles for five air showers emerging from the ground at depths of 10000, 25000, 50000, 
75000, and 100000 g/cm 2 (dashed lines). 
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Fig. 2. — The calculated HiRes-II electron-neutrino aperture. 
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Fig. 3. — The H iRes-H neutrin o flux limit, black boxes: v e limit (this work), cross-hatched 
boxes: v T limit (jMartena I20071 ). open boxes: v e and u T combined flux limi t. Dotted line : 
cosmogenic per flavor neutrino flux limit from fits to HiRes cosmic-ray data (jBrusovall20071 ). 
Dashed line: cosmog enic per flavor neutrin o flux limit derived from fits to existing cosmic- 
and gamma-ray data (jSemikoz fc Siglll2004l ). Dot-das hed line: cosmogenic p er flavor neutrino 



flux from fits to HiRes and AG AS A cosmic-ray data (|Seckel fc Stanev 



calc ulated neutrino flux limits fro m the Fly's Eve (jBaltrusaitis et al. 



2005). Also shown are 



1984. 



1985), A NITA- 



lite JOorham et allbooil. R ICE kravchenko et al.lboOfil ). AGASA jchikawa et all boOll ) 



and Auger ([Abraham et al. 1120081 ) experiments. 
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